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Abstract 

We present the results of a lattice QCD calculation of the pseudoscalar meson 
decay constants f n , fx, /d and fo 3 , performed with Nf = 2 dynamical fermions. 
The simulation is carried out with the tree-level improved Symanzik gauge action 
and with the twisted mass fermionic action at maximal twist. We have considered 
for the final analysis three values of the lattice spacing, a ~ 0.10 fm, 0.09 fm and 0.07 
fm, with pion masses down to ~ 270 MeV. Our results for the light meson decay 
constants are fx = 158.1(2.4) MeV and fx/In = 1.210(18). From the latter ratio, 
by using the experimental determination of T(K — > ^^(7)) /T(ir — > ^^(7)) and the 
average value of \V UC [\ from nuclear beta decays, we obtain \V US \ = 0.2222(34), in good 
agreement with the determination from semileptonic K13 decays and the unitarity 
constraint. For the D and D s meson decay constants we obtain fjj = 197(9) MeV, 
fo a = 244(8) MeV and /d s //d = 1-24(3). Our result for fjj is in good agreement 
with the CLEO experimental measurement. For fo s our determination is smaller 
than the PDG 2008 experimental average but in agreement with a recent improved 
measurement by CLEO at the 1.4 a level. 

1 Introduction 

An accurate lattice determination of the pseudoscalar decay constants of kaon and D- 
mesons is an important task. On the one hand, the ratio fx / fn can be used together with 
the experimental determination of T(K — > jiu^ij)) /T(ir — > /iz/^7)) and the average value 
of \V u d\ from nuclear beta decays to achieve a precise determination of the CKM matrix 
element |V^ S | [1] and to test the CKM first raw unitarity relation. On the other hand, the 
pseudoscalar decay constants fr> and fr> s have been recently measured at CLEO, BaBar 
and Belle [21 E], asking for accurate lattice determinations to be compared to. 

In this paper, we present a lattice QCD calculation of the pseudoscalar meson decay 
constants f n , fx, fo and fc s - With respect to our previous study of the pion and kaon 
decay constants [4], here we have analysed data at three values of the lattice spacing, 
a ~ 0.10 fm, 0.09 fm, 0.07 fm (corresponding to j3 = 3.8,3.9,4.05), and performed a chiral 
extrapolation taking lattice artefacts into account. Estimating the lattice artefacts turns 
out to be crucial for an accurate determination of fo and fu s since cutoff effects induced 
by the charm mass, which are parametrically of 0(a 2 m 2 c ), are not small in our simulation, 
~ 5 4-10%. Both SU(2) and SU(3) chiral perturbation theory (ChPT) has been considered 
for the chiral extrapolation, whereas only the latter was considered in [4J. With respect to 
ref. jl], we have also added ensembles with a lighter quark mass (m^ ~ 270 MeV) and a 
larger volume (L ~ 2.7 fm), both at the value of a ~ 0.09 fm. 

The calculation is based on the gauge field configurations generated by the European 
Twisted Mass Collaboration (ETMC) with the tree-level improved Symanzik gauge ac- 
tion [5] and the twisted mass action [6] at maximal twist, discussed in detail in refs. [7]-|10j. 
We simulated Nf = 2 dynamical quarks, taken to be degenerate in mass, whose masses are 
eventually extrapolated to the physical isospin averaged mass of the up and down quarks. 
The strange and charm quarks are quenched in the present calculation. 

The use of the twisted mass fermions turns out to be beneficial, since the pseudoscalar 
meson masses and decay constants, which represent the basic ingredients of the calculation, 
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are automatically improved at 0{a) [TT] (see also [12]), and the determination of the 
pseudoscalar decay constants does not require the introduction of any renormalization 
constant. Both these features allow to significantly improve the accuracy of the calculation. 
It has also been shown that, for twisted mass fermions at maximal twist, the so called 
KLM factor [13J, which relates the lattice quark propagator at zero momentum to the 
continuum one, is equal to one at tree level, to all order in am q [H]. This is beneficial 
in reducing discretization effects particularly for heavy quark masses, as the charm quark 
mass considered in this study. 

As discussed in refs. (U [151 IB], we implement non-degenerate valence quarks in the 
twisted mass formulation by formally introducing a twisted doublet for each non-degenerate 
quark flavour. In the present analysis we thus introduce in the valence sector three twisted 
doublets, (u, d), (s, s') and (c, c'), with masses /xj, /i s and /i c respectively. Within each dou- 
blet, the two valence quarks are regularized in the physical basis with Wilson parameters 
of opposite values (r = — r' = 1). We simulate mesons composed of quarks with opposite 
Wilson parameters so that the squared meson mass m 2 PS differs from its continuum coun- 
terpart only by terms of 0(a 2 n q ) and 0{a i ) ) whereas fpg differs from its continuum limit 
by terms of 0(a 2 ) [T7J [18]. This implies that at 0(a 2 ) the cutoff effects on m 2 PS and fps 
are as in a chiral invariant lattice formulation. In our calculation large artefacts, like those 
affecting the neutral pion mass in the twisted mass formulation of lattice QCD, seem not 
to be present. 

In the present analysis, we study the pseudoscalar decay constants as a function of 
the meson masses, whereas in our previous work [I] we relied on their dependence on 
the quark masses. When data at different values of the lattice spacing are involved, the 
study in terms of meson masses is simpler, since it does not require the introduction of 
the quark mass renormalization constant [Z m = Zp 1 ) to convert at each lattice spacing 
from the bare to the renormalized (cutoff independent) quark mass. The dependence of 
the decay constants on the meson masses is studied together with the dependence on the 
lattice spacing, through a combined fit where terms of 0(a 2 ), coming from the Symanzik 
expansion of the lattice theory, are added to the functional forms predicted by ChPT. 
In this way, the continuum and chiral extrapolations of the lattice results are performed 
simultaneously. In alternative to this combined analysis, a different approach could be 
adopted. It consists of first extrapolating data at fixed meson mass values to the continuum, 
and then extrapolating the obtained continuum results to the physical point. With our 
simulation setup, however, this procedure turns out to be unsafe, since for some values of 
the meson masses we have data at only two values of the lattice spacing. Such a procedure 
could become feasible when data at a smaller value of the lattice spacing are available. 
Corresponding simulations with a ~ 0.05 fm are currently performed by ETMC 

In order to perform the extrapolation to the physical masses we have used ChPT for 
the light mesons and Heavy Meson ChPT (HMChPT) [19] for the D sector. For the kaon 
and D s mesons we have considered both SU(2)- and SU(3)-ChPT. In the SU(2) case one 
treats the u/d quarks as light, while the strange quark is not required to satisfy chiral 
symmetry. The short interpolation to the physical strange quark, which is required in our 
analysis, is then performed linearly [20J. This is justified, since our simulated values of the 
strange quark mass are around the physical mass. For comparison, we have also considered 
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chiral extrapolations based on SU(3)-ChPT and SU(3)-HMChPT. We find that the SU(2) 
effective theory, which is less predictive than SU(3), provides however a better description 
of the lattice data for the decay constants up to the region of the strange quark mass. 

Our final results for the kaon and D-mesons decay constants are given in the abstract 
and in eqs. ( ITT]) and ( 1231) . Our determination of fx/ fn leads to a determination of \V US \ 
that is in good agreement with the value obtained from semileptonic kaon decays, though 
with a larger error, as well as with the first row unitarity constraint of the CKM matrix. Of 
relevant phenomenological interest is also our result for fjj a , which is about 2.3 a lower than 
the experimental average quoted by the PDG [2] but in agreement with other unquenched 
lattice determinations and with a recent improved measurement by CLEO [3] at the 1.4 a 
level. The value indicated by the new CLEO measurement weakens the tension between 
experimental and lattice results for fo s , which suggested explanations in terms of new 
physics effects [2T| 122]. 

The plan of this paper is as follows. In section 2 we provide the details of the lattice 
simulations used for the present study and discuss the determination of the pseudoscalar 
meson masses and decay constants. The combined chiral and continuum extrapolation fits 
of the light and D-mesons decay constants are discussed in sections 3 and 4, respectively. 
There, we also provide our final results for the decay constants, discussing in particu- 
lar the evaluation of the systematic uncertainties and the comparison with other lattice 
determinations and with recent experimental measurements. 

2 Simulation details 

Details of the ensembles of gauge configurations used in the present analysis and the values 
of the simulated valence quark masses are collected in Tables [TJ and [2J respectively. 

Measurements are performed over independent gauge configurations that are separated 
by 20 trajectories in the case of the ensembles at (3 = 3.8 and 3.9, and by 40 trajectories 
in the case of the ensembles at (3 = 4.05. The trajectory length is equal to unity for the 
ensembles A 2 , A s and B 7 and to 1/2 for the ensembles B1-B4, B 6 and C1-C3. Among the 
available ETMC ensembles we have excluded from this study those corresponding to pion 
masses larger than 500 MeV (ensembles A4, B 5 and C4 of ref. [9]). 

Our strategy and the conditions to tune to maximal twist have been discussed in refs. [7]- 
[9]. Whereas at (3 = 3.9 and (3 = 4.05 these conditions are accurately fulfilled, at (3 = 3.8 
the situation is different. Due to large fluctuations and long autocorrelations for the PCAC 
mass appearing at the smallest value of the twisted mass parameter at (3 = 3.8, we cannot 
be confident that the maximal twist condition is realized with the same accuracy as at 
(3 = 3.9 and [3 = 4.05. In the present study, in order to check for effects of such a possible 
mismatch, besides not considering the ensemble with the lightest quark mass at (3 = 3.8 
(ensemble Ai of ref. [9]), we have also performed an analysis with and without taking the 
whole set of data at (3 = 3.8 into account. As we will demonstrate below, fully consistent 
results are obtained. This finding is also supported by the results of a theoretical analysis 
of the effects of being out of maximal twist, which shows that these systematic effects on 
the pseudoscalar decay constants analyzed here at (3 = 3.8 are small compared to statistical 
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/3 

P 


a [rmj 


V/a 




m n [MeVJ 


m n L 


AT 

N cfg 






A 2 


3.8 


0.10 


24 3 x 48 


0.0080 


410 


5.0 


240 


[10,23] 


[14, 23] 


A 3 








0.0110 


480 


5.8 


240 






B x 


3.9 


0.085 


24 3 x 48 


0.0040 


315 


3.3 


480 


[12,23] 


[16, 23] 


B 2 








0.0064 


390 


4.0 


240 






B 3 








0.0085 


450 


4.7 


240 






B 4 








0.0100 


490 


5.0 


240 






B 7 


3.9 


0.085 


32 3 x 64 


0.0030 


270 


3.7 


240 


[12,31] 


[16,31] 


B 6 








0.0040 


310 


4.3 


240 






Ci 


4.05 


0.065 


32 3 x 64 


0.0030 


310 


3.3 


144 


[16,31] 


[21,31] 


c 2 








0.0060 


430 


4.6 


128 






c 3 








0.0080 


500 


5.3 


128 







Table 1: Details of the ensembles of gauge configurations used in the present study: value 
of the gauge coupling (3; approximate value of the lattice spacing a; lattice size V = L 3 xT 
in lattice units; bare sea quark mass in lattice units; approximate value of the pion mass; 
approximate value of the product m n L; number of independent configurations N c f g . We 
also provide the htting time interval in lattice units chosen for the two-point pseudoscalar 
correlators in the pion and kaon sectors, (At)^^, and in the D -meson sectors, (At)n,D s - 





A 2 -A 3 


B\ — B 4 


B 7 


B e 


Ct-C 3 


a/j,i 


0.0080, 0.0110 


0.0040, 0.0064, 0.0085, 
0.0100 


0.0030 


0.0040 


0.0030, 0.0060, 
0.0080 




0.020, 0.025, 0.030, 
0.036 


0.022, 0.027, 0.032 


0.022, 
0.027 


0.022, 
0.027 


0.015, 0.018 
0.022, 0.026 


a/j, c 


0.270, 0.310, 0.355, 
0.435, 0.520 


0.250, 0.320, 0.390, 
0.460 


0.250, 
0.320 


0.250, 
0.320 


0.200, 0.230, 
0.260 0.315 



Table 2: Values of simulated valence quark masses in lattice units for each configuration 
ensemble in the light, strange and charm sectors. 



and other systematic uncertainties on the same datajj 

With respect to our previous determination of /„- and fx jl], the new ensembles used in 
the present analysis are those with (3 = 3.8 (A 2 -A 3 ) and (3 = 4.05 (Ci-C 3 ) and the ensembles 
Bq-Bj at (3 = 3.9. The ensemble B-j has the lightest simulated mass fii ~ 0.15 mP s hys ', where 
m phys. j g phy S i ca i strange quark mass. The ensembles B\ and Bq have the same value 
of (3 and sea quark mass but different volumes, L ~ 2.0 fm and L ~ 2.7 fm, thus allowing 
a study of finite size effects. 

In order to investigate the properties of the ir, K, D and D s mesons, we simulate 
the sea and valence light (u/d) quark mass in the range 0.15 m^ hys ' ^ fii ^ 0.5mf hys ', 
the valence strange quark mass is within 0.9mf !/s ' ^ fi s ^ 1.5 m p J iys ', and the valence 
charm quark mass within 0.8™^' ^ /i c ~ 1.5m phys ', m p c hys ' being the physical charm 
mass. The values of the valence quark masses simulated for each configuration ensemble 

1 For the basic ideas of this analysis in the unitary case see ref. [25] . 
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are collected in Table [2j These values can be converted into the corresponding values of 
the renormalized quark masses (e.g. in the MS scheme) using the available determinations 
of renormalization constants given in refs. [2H [25] . 

As already noted, with twisted mass fermions at maximal twist, the determination of 
the pseudoscalar decay constants, besides being automatically improved at 0(a), does not 
require the introduction of any renormalization constants. For a pseudoscalar meson of 
mass mps-, composed of valence quarks with masses /i^ and //^, the decay constant fp$ 
is obtained as 



where P = qij^- The meson mass mps and the matrix element |(0|P|PS)| entering 
eq. ([Q) have been extracted from a single state fit of the two-point pseudoscalar correlation 
function within the time intervals collected in Table [TJ The replacement of mps with 
sinh mps in the lattice definition of the decay constant helps in reducing discretization 
errors for heavy meson masses. Combined with the observation that the tree-level KLM 
factor for the quark field is equal to one for twisted mass fermions at maximal twist [14] , 
this replacement allows to remove at tree level all 0((am c ) n ) terms in the determination 
of the P( s )-meson decay constant. 

The statistical accuracy of the meson correlators is improved by using the so-called "one- 
end" stochastic method, implemented in ref. [26] (see also pj), which includes all spatial 
sources at a single timeslice. Statistical errors on the meson masses and decay constants 
are evaluated using the jackknife procedure, with 16 jackknife bins for each configuration 
ensemble. Statistical errors on the fit results which are based on data obtained from 
independent ensembles of gauge configurations are evaluated using a bootstrap procedure, 
with 100 bootstrap samples. In order to illustrate the quality of the data, we show in 
fig. [T] the effective masses of pseudoscalar mesons, as a function of the time, for four 
choices of quark mass combinations, representing the pion, the kaon, the D and D s mesons, 
respectively. The pseudoscalar masses shown in fig. [1] are extracted from the two point 
correlator of the charged pseudoscalar density, together with the matrix element (0\P\PS). 
Both quantities enter eq. ([1]) for fp$. Each plot shows the effective mass obtained from the 
ensembles B\ and P 6 , i.e. at /3 = 3.9, with a\i sea = 0.0040 and with two different lattice 
sizes, namely 24 3 x 48 and 32 3 x 64 (see Sec. 3.1 for a discussion of finite size effects in our 
analysis) . 

3 The pion and kaon decay constants 

3.1 Combined chiral and continuum extrapolation 

A good convergence of SU(3)-ChPT is in general not guaranteed in the kaon sector. As 
recently pointed out in [20J (see also [27] for a detailed review on this subject), a safer 
approach consists in avoiding the chiral expansion in terms of the strange quark mass and 
applying therefore SU(2)-ChPT. The use of SU(2)-ChPT is well motivated in our analysis, 
since we simulated fx s in the range of the physical strange quark mass, thus having small 
values of \i s (see Table At next-to-leading order (NLO), the SU(2)-ChPT prediction 
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Figure 1: Effective pseudoscalar meson masses rrfpg (a fi sea , a fx]^, a fi^) as a function of 

(1 2) 

time, in lattice units, with /j, seci and fi^ denoting genetically the sea and valence quark 
masses respectively. For illustrative purposes the following choices of quark mass com- 
binations are displayed: m P5 (0.0040, 0.0040, 0.0040) (pion), m P5 (0.0040, 0.0040, 0.0220) 
(kaon), m P5 (0.0040, 0.0040, 0.2500) (D-meson), m P5 (0.0040, 0.0220, 0.2500) (D s -meson). 
In each plot we compare the effective masses as obtained from the two ensembles B\ and 
Bq, which correspond to different lattice sizes. Dashed and solid lines represent the 1-a 
ranges of the corresponding masses as obtained from the fit of the two-point correlation 
functions. 



for the pion decay constant is well known [28J, 

fps(jH,fii,fJti) = f ■ (1 - 2&]n& + &&) , (2) 
and the corresponding expression for the kaon decay constant reads [20J 



fpsfa, IH, I*.) = f (K) ■ (l - + . (3) 

We are using the notation (fi sea , //^, fi^) for the quark mass content of the corresponding 
meson, and the variables £'s in eqs. (j2j) and ([3]) are expressed in our analysis as a function 
of meson masses J! 



2 We use the normalization in which f„ = 130.7 MeV. 
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The leading contribution in eq. ([2]) is represented by the low energy constant (LEC) /, 
which provides the value of the pion decay constant in the chiral limit, whereas the coeffi- 
cient b is related to the LEC I4 of the NLO chiral Lagrangian. Notice that, in the Nf = 2 
theory we are simulating, i.e. with a quenched strange quark, these constants are indepen- 
dent of the value of the strange quark mass. In this theory there is actually no distinction 
between SU(2) and SU(3) ChPT expansions for pion observables. The LECs f^ K > and 
entering the SU(2) formula d2J) for the kaon decay constant, instead, are functions of the 
(valence) strange quark mass. 

In order to perform a combined fit of the data for the pseudoscalar decay constants at 
the three values of the lattice spacing, we rely on the Symanzik expansion of the lattice 
regularized theory and introduce in the fitting formulae discretization terms of 0(a 2 ) and 
0(a 2 fig). Discretization effects of C(a 2 /i/), i.e. proportional to the light quark mass, 
represent very small contributions that turn out to be invisible in the fit. Moreover, since 
the simulated /i s masses are all close to the physical strange quark mass, we can safely 
linearize the strange mass dependence of the LECs f( K > and U- K ' in eq. (jSJ) around m p s hys \ 
We thus write the SU(2)-ChPT fitting formulae for the pion and kaon decay constants as 

f P3 (jii,Hi,fii) = f-(l-2£ u )n£u + b£u + A^), (5) 



, 



1 - -i u ln& + + C } U + (A a + A as £ ss )^ 
4 r 



(6) 



The fit is performed in units of the Sommer parameter tq [29]. The values of r^ja at 
the three lattice spacings have been extracted in ref. [9] from the analysis of the static 
potential, obtaining 

r /a = {4.46(3), 5.22(2), 6.61(3)} (7) 

at (3 — {3.8, 3.9, 4.05} respectively. The physical value of tq in the continuum limit is 
determined in our analysis by combining the determination of tq ■ f n with the experimental 
value of the pion decay constant. This procedure, combined with eq. (j7j), corresponds to 
fixing the lattice scale using f n as physical input. 

An important ingredient in the analysis is the study of finite size effects (FSE). With 
our simulation setup, the largest FSE are expected in the data of the ensembles B\ and 
Ci, for which m n L ~ 3.3 (see Table [1]). A quantitative estimate of these effects can be 
obtained from the comparison of the data of the ensembles B\ and B 6 , that only differ in 
lattice size. This comparison provides consistent results with the FSE predicted at NLO 
by SU(2)-ChPT [301 EI], which are expressed for the pion and kaon decay constants by 

fps(fii, Hi, M L) = fpsi.Hu m, Hi) ■ [1 -2£u<7i(L,£h)] , 

3 

fps(Vh V-h Vs] L) = fpsim, Hi, Hs) • 1 - 7 iu 9i(L, £u) 

where the function gi is defined for instance in ref. @]@ This correction, which on the 
ensembles B\ and C\ corresponds to about 2.5% for f n and 0.9% for fx, has been included 

3 Note that the finite size corrections in eq. ([5]) are obtained from the loop contribution of the infinite 
volume ChPT predictions in eqs. @ and by replacing ln£ with the function <?i(L,£). 
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ro ■ / 


r o ■ Jo 


ro ■ /if) 


b 


°0 


b (K) 


A 


A a 


A a s 


0.271(6) 
0.274(9) 


0.305(6) 
0.312(9) 


0.18(1) 
0.18(2) 


-0.25(13) 
-0.33(15) 


0.5(1) 
0.4(2) 


-1.9(4) 
-1.7(5) 


0.7(5) 
0.5(7) 


0.7(4) 
0.3(7) 


2.1(6) 
3(1) 



Table 3: Values of the SU(2) fit parameters of eqs. and (0), as obtained by including 
(first row) or excluding (second row) the data at (3 = 3.8. Quoted errors are statistical 
plus fitting errors. 





r [GeV" 1 ] 


fx [MeV] 


JkI fn 


X'fdof 


SU(2)-ChPT 


incl. 13 = 3.8 
excl. /? = 3.8 


2.22(5) 
2.25(7) 


158.1(8) 
158.9(1.4) 


1.210(6) 
1.216(11) 


{11/8; 40/30} 
{7/6; 35/22} 


SU(3)-ChPT 


incl. (3 = 3.8 
excl. (3 = 3.8 


2.23(5) 
2.28(7) 


158.2(6) 
158.0(0.8) 


1.210(5) 
1.209(6) 


61/42 
54/32 



Table 4: Values of tq, Jk and fx / fir as obtained from SU(2)- and SU(3)-ChPT by including 
or excluding the data at (3 = 3.8. For each fit, the chi-squared per degree of freedom, 
X 2 /dof, is also given in the last column. For fits based on SU(2)-ChPT, the two values of 
X 2 /dof correspond to the fit of f n and fx respectively. Quoted errors are statistical plus 
fitting errors. 



in our fit. For a more detailed discussion of FSE in the pion decay constant see ref. [8]. In 
our data for the kaon decay constant, instead, the differences between the ensembles 73 1 and 
B e are small and at the level of the statistical errors. We obtain Af K = f K =24 /f K =32 ~ 1 = 
0.005(4), compatible with zero. Also the NLO partially quenched ChPT prediction for 
A/ x [31] is similarly small (—0.006), at the level of our statistical error. 

The values of the fit parameters and the results for r , f K and fx/ fn are collected in 
Tables [3] and HI respectively. The extrapolation to the physical up/down quark mass and 
the interpolation to the physical strange mass has been performed by inserting in eqs. ([5]) 
and © tu = (mP hys -) 2 /(Anf) 2 and £ ss = (2 {m p ^ ys -f - (mf ys -) 2 )/^7rf) 2 , where m^ ys - and 
m phys. are experimental pion and kaon masses. In both Tables, we show the results of 
our fits when we take the data at [3 = 3.8 into account and when we leave them out. As 
can be seen, the values of the fit parameters, as well as those of the decay constants, are 
found to be well consistent in the two cases. In the following, therefore, we will consider 
for fx and fx/ f-w only the predictions obtained by including the (3 = 3.8 data. From 
the results given in Tables [3] and HI one can also derive our prediction for the pion decay 
constant in the chiral limit, /, and the LEC I4 = b/2 + 2 ln(47r//m 7r +). We obtain the 
values / = 121.7(l)MeV and U = 4.66(6), which are in good agreement with the results 
of the scaling analysis performed by our Collaboration in [10J, / = 121.66(7) (26)MeV and 
I4 = 4.59(4)(13). The quality of the fit for the combined chiral and continuum extrapolation 
of the pion and kaon decay constant is illustrated in fig. [2J 

As an alternative to the SU(2)-ChPT approach, we have also considered the expansion 
valid for a small strange quark mass, fitting both the pion and the kaon decay constants 
using SU(3)-ChPT. The relevant expression, valid for the partially quenched Nf = 2 theory 
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Figure 2: Lattice results for r f n = r Q f PS (ni, n h ni) and r f K = r f PS (fi h y, h fi s ) as a 
function of the pion mass square {r-Qm^) 2 = {rompsinu^h^i)) 2 - F° r the kaon, we dis- 
play data with fi s fixed to the simulated mass that corresponds to a reference meson mass 
TQmpsini-, /i s , — 1.63. The SU(2)- (SU(3)-) ChPT extrapolation to the physical pion 
mass is represented at fixed lattice spacing by the dashed (dotted ) curves, and in the con- 
tinuum limit by the solid (dashed-dotted) curve. Our results for the physical values of the 
pion and kaon decay constants, obtained from SU(2)-ChPT, are illustrated by diamonds in 
the plot. In the kaon case an interpolation to the physical strange quark mass is performed. 

at NLO, is [32] 



fps(f*h to, Us) = f- 



2 £h 

1 - -in In & - — ln£ ss - £, a ln& s + b u £ u + b ss £ ss + (B a + B as £ s 



r 2 



(9) 

where, as in the SU(2) case, we have also included in the fit discretization terms of 0(a 2 ) 
and C(a 2 /i s ) as well as finite size corrections [3T] . 

The results obtained from the SU(3)-ChPT analysis are given in Table H] and shown in 
fig. [21 and are found in very good agreement with those obtained from the SU(2) fit. A 
more careful analysis suggests, however, that the SU(3)-ChPT fit is less robust than the one 
based on SU(2). In SU(2)-ChPT, at NLO, one obtains a good fit of the data by expressing 
the chiral formulae (jSJ) and ([6]) either in terms of meson masses, as performed here with 
£ij defined as in eq. PJ, or in terms of quark masses, i.e. with £y = B (/ij + /ij)/(47r/) 2 , 
where B is the LEC entering at LO in the chiral Lagrangian. In SU(3)-ChPT, instead, 
the NLO formula expressed in terms of meson masses provides a good description of the 
lattice data, but fits performed in terms of quark masses require in the kaon sector the 
inclusion of NNLO terms, as we already found in ^j. This means that the replacement 
of quark with meson masses effectively resums higher order chiral contributions, actually 
improving the fit based on NLO SU(3)-ChPT of the pseudoscalar decay constant beyond 
mps ~ 450 Me V. A similar result was found in ref. |33j . 
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3.2 Results for fx and JkI 

As seen in the previous section, lattice data in the kaon sector could be analysed by means 
of either SU(2)- or SU(3)-ChPT, leading to almost identical results, see TableSJ The results 
we quote for f K and fx/fn are those obtained from SU(2). The errors quoted in Table H] 
are statistical plus fitting errors from the combined chiral and continuum extrapolation. 
We now discuss how we evaluate other systematic errors. 

Since we have simulated at three values of the lattice spacing and on our coarsest 
lattice (/3 = 3.8) we have data for only two values of the light quark mass, we include in 
our final results a systematic uncertainty due to residual discretization effects. The leading 
discretization errors in our determination of the light meson decay constants are expected 
to be of 0(a 2 Aq CD ). This naive expectation is roughly confirmed by the results of our fit. 
On our finest lattice, for instance, with (3 = 4.05 and a ~ 0.07 fm, one has a 2 A.q CD ~ 1-^2% 
and we find that the difference between the values taken by the kaon decay constant on 
this lattice and its estimate in the continuum limit is approximately 2.6% (this difference 
turns out to be slightly larger, about 2.8%, at the reference mass r m PS (fi[, fi s , [i s ) = 1.63 
for which results are displayed in fig. [2]). We conservatively assume an error of 50% in the 
continuum extrapolation starting from the data on our finest lattice, thus adding to our 
final results for fx and fx:/ fit a relative systematic error of 1.3% (half of the difference 
between the values at /3 = 4.05 and the continuum estimates). 

In the present analysis, FSE corrections have been implemented by using the predictions 
of NLO ChPT, as discussed in the previous section. Besides the direct comparison of this 
theoretical estimate with the data available on the two lattices B\ and B e , where m n L 
varies from 3.3 to 4.3, an additional indication that these corrections are under control 
is provided by the compatibility between the results for f v determined here, by treating 
the FSE with NLO ChPT, and those obtained in [S] by using the resummed formulae of 
ref. [3l]. For the kaon decay constant, FSE are found at the level of the statistical errors 
at most. A fit performed without including this correction provides a result for the kaon 
decay constant which is lower by about 0.7%. We conservatively include this difference in 
the systematic error of fx as an estimate of the uncertainty due to FSE. 

The only uncertainty which cannot be reliably estimated within our Nf = 2 simulation, 
is the error due to the quenching of the strange quark. The good agreement observed among 
the recent Nf = 2 and Nf = 2 + 1 lattice determinations of fx and fx/ fn (see fig. [3]) 
suggests, however, that such an effect is smaller than the other systematic uncertainties 
estimated above. The ETM Collaboration is planning to investigate directly the effect 
of the quenching of both the strange and charm quarks through simulations with Nf = 
2 + 1 + 1, which are currently in progress [35]. For a more extensive discussion of the 
various sources of lattice systematic uncertainties see ref. [36] . 

Our final results for the kaon decay constant and the ratio fx / fn are then 

fx = 158.1(0.8)(2.0)(1.1) MeV , f K /f v = 1.210(6) (15) (9) , (10) 

where the first error comes from statistics and chiral extrapolation, the second from the 
estimate of residual discretization effects and the third from the uncertainty on FSE cor- 
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Figure 3: Lattice QCD determinations of the ratio fx/fn obtained from simulations with Nf = 
2 \37\ and Nf = 2 + 1 {23 J37 ) /-J43 y dynamical quarks. A star in the legend denotes preliminary 
results. The results are also compared with the experimental average of fx / fn obtained by using 
for V us the determination from decays \44j . 



rections. By combining the errors in quadrature, we finally obtain 

f K = 158.1(2.4) MeV , f K /f w = 1.210(18) . (11) 

Our result for the ratio fx/fn is compared in fig. [3] with other unquenched lattice 
determinations, performed with either Nf = 2 or Nf = 2 + 1 dynamical quarks, as well 
as with the experimental average of fx/ fir obtained by using for V us the determination 
from K.£ 3 decays [Hj. Our result turns out to be in very good agreement with the latter 
determination, as well as with most of the Nf = 2 and Nf = 2 + 1 lattice results. 

Alternatively, our result for fx / fn can be combined with the experimental measurement 
of T(K — > /iP At (7))/r(7r — ► /iz/ Ai (7)) [44] to get a determination of the ratio |V^ s |/|V^d| PP- 
We obtain 

|K.|/|Kd| = 0.2281(5)(35), (12) 

where the first error is the experimental one and the second is the theoretical error com- 
ing from the uncertainty on fx/ fn- Eq. (I12j) . combined with the determination \V u d\ = 
0.97425(22) [45J from nuclear beta decays, yields the estimate 

| V us \ = 0.2222(5) (34), (13) 

in good agreement, though with a larger error, with the value extracted from decays, 
\V US \ = 0.2246(12) [S]. Eq. flEJJ and the value of \V ud \ quoted above leads to 

\V ud \ 2 + \V US \ 2 + \V ub \ 2 - 1 = -1.5(1.6) ■ 10- 3 . (14) 

in good agreement with the unitarity constraint of the CKM matrix. 



11 



4 The D and D s decay constants 



4.1 Combined chiral and continuum extrapolation 

In order to determine the D and D s meson decay constants we essentially proceed as in the 
kaon sector. We analyse simultaneously data at the three values of the lattice spacing and 
perform for the pseudoscalar decay constants combined fits of the meson mass dependence 
and of discretization effects. The simulated values \x c of the charm quark mass are close to 
the physical charm quark mass (0.8 m v ^ ys ' ^ /i c ^ 1.5m^ ys '), so that the interpolation to 
the physical value is short and smooth. From the comparison of the data of the ensembles 
B\ and B§ we also find that FSE are negligible for the D s decay constants, and they are at 
the level of the statistical error or smaller for fp. On the other hand, discretization errors 
induced by the charm quark mass have to be taken carefully into account in the fit, being 
parametrically of 0(a 2 i.e. approximately 5 -r- 10% in our simulation. 

The functional forms describing the mass dependence of the decay constants assumed 
to fit the data in the D and D s sectors are those predicted by HMChPT [TH] . We consider 
the SU(2) version of the theory, as in the case of the kaon sector, where the strange quark 
is not required to satisfy chiral symmetry, but it is considered heavy enough to justify an 
expansion in powers of fii/ '/j, s . For comparison, we have also investigated the predictions 
of SU(3)-HMChPT where, instead, the strange quark is required to satisfy the same chiral 
symmetry of the light up and down quarks. As we will see below, in the D-meson sector 
the SU(2)-HMChPT approach turns out to work significantly better than the one based 
on SU(3). 

Within the SU(2)-HMChPT analysis, we extract fjj and fo s by considering two different 
procedures. In the first one we fit the two following combinations of meson masses and 
decay constants: 

f Ds ^, R= f ° sV ™^ , (Fit I). (15) 

JD\J™D 

We find, in particular, that the advantage of introducing the ratio R is that discretiza- 
tion effects largely cancel in the ratio. In the second approach we consider the previous 
quantities divided by the light decay constants, i.e. we fit the ratios 

Rl = fD.v**: ^W^r x _^. (FitII) (16) 

Ik Jk jDy/mB 

Here, the advantage of the ratio i?2 is that it exhibits a quite smooth chiral behaviour. The 
comparison of the results obtained for fa, fr> s and /dJ Id in the two cases, Fit I and Fit 
II, will provide an estimate of the systematic uncertainty due to the chiral extrapolation. 

In all quantities entering eqs. (|T5|) and ([IB]) , the D-mesons decay constants are mul- 
tiplied by the square roots of the corresponding meson masses, in order to reconstruct 
the observables that remain finite in the infinite mass limit. The Heavy Quark Effec- 
tive Theory predicts in fact for a Heavy (if)-light(Z) meson an expansion of the form 
fm-s/TriHi = A + B/ttihi + C(V m H/)) U P to small radiative corrections. Though the 
heavy quark expansion is known to be slowly convergent in the charm mass region, in 
our analysis we can safely assume such a behaviour for the D mesons since only a short 
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interpolation of the lattice data to the physical charm quark mass is needed. Moreover, 
since the contribution of the sub-leading 1/mm correction in this interpolation is small, 
we can safely account for the dependence on the light meson masses only in the leading 
term, by using the prediction of HMChPT. 

We obtain the SU(2)-HMChPT functional forms for the quantities in eqs. (1151) and ( [161) 
by expanding the corresponding SU(3)-HMChPT predictions in powers of /x;//i s , and re- 
absorbing the strange quark mass dependence in the SU(2) LECs0 The SU(3)-HMChPT 
prediction for f Dg ^/m Da , valid in the partially quenched Nf = 2 theory [19], is given by 



3/2 
r 



1 + 3^ c 2 



til ~ 2£ 



In 



lnf M + C 2 £ii + C£ 8 



+ 



Ca 







(17) 

where the parameter g c is related to the gD*D-n coupling by go*D-K = {2^/mr>mr>* / f^) g c . 
Since eq. (I17p does not contain logarithms of the pion mass (i.e. ln^), one finds that 
its expansion in powers of \i a — Cu/Css leads to an SU(2) chiral expression for the D s 
decay constant which is free of chiral logarithms at NLO. As in the light meson case, we 
also include in the fitting formula discretization terms of 0(a 2 ), in order to take simulta- 
neously into account the lattice artefacts. We observe, in this respect, that the Symanzik 
expansion of fjj and fo s contains at 0(a 2 ) discretization terms depending on the charm 
quark mass either linearly or quadratically, with the leading contribution expected from 
terms of 0(a 2 fi 2 ). The limited set of data available in our analysis, however, does not 
allow us to fit both these dependencies separately. We parameterize these discretization 
effects in terms of the meson masses, and we thus introduce in the fitting formula only a 
term proportional to a 2 m 2 D . We have also tried an alternative fit where the charm mass 
dependent discretization term is taken to be proportional to a 2 m Ds instead of a 2 m 2 Ds , and 
we obtain completely consistent results. Thus, we use as our fitting formula for fo s ^/wifl, 
the expression 



3/2 



1 + D 2 Zu + (D a + D a £ ss ) — + D ah a 2 m 2 Ds 



+ 



D, 



5/2 ' 

r m D 3 



where the coefficients Di and D 2 , which depend on the strange quark mass, are expressed 
in the fit as linear functions of this mass: 



Di — Difi + Di !m £ s 



(19) 



with z = l,2. 

The SU(2)-HMChPT prediction for the ratio R, defined in eq. (fTBT) . is straightforwardly 
obtained by dividing eq. ( fl8l) by the HMChPT expression for /ci/mo, which is provided 
by the SU(3)-HMChPT formula of eq. ( JT7I) with fi s = n\. Thus, we assume for the ratio R 
the expression 

,2 



R = D[ 



1 + - (l + 3g 2 ) £« In £« + D'^i + (D' a + D' as U ^ + D' a 



2 2 
ah a m Ds 



+ 



D' 



r m Ds 



(20) 



4 The same procedure allows to obtain the SU(2)-ChPT expression © of the kaon decay constant fx 
from the SU(3) prediction of eq. ©• 
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^1,0 




^2,0 




n„ 

^3 


n 


n 


Uah 


1.62(9) 
1.54(9) 


0.78(7) 
0.9(1) 


0.4(2) 
0.7(2) 


-0.7(5) 
-1.3(5) 


-2.9(3) 
-2.6(3) 


-0.2(4) 
-0.6(7) 


-2(1) 
-3(2) 


0.13(3) 

0.16(2) 




rv 


^2,0 




^3 




D'as 


D'ah 


1.0(1) 
1.0(1) 


1.1(2) 
1.2(2) 


1.9(4) 
2.1(6) 


-2.1(9) 
-3(1) 


0.6(4) 
0.5(6) 


-1.4(9) 
-1(1) 


0(1) 

1(2) 


0.07(5) 
0.07(7) 




A,o 


P\.m 


-^2,0 


P2,m 


P 3 


Pa 


p 

- 1 as 


Pah 


4.9(2) 
4.7(2) 


0.7(2) 
0.8(3) 


0.7(2) 
0.9(2) 


0.1(5) 
-0.5(5) 


-7.2(8) 
-6.8(7) 


-0.1(4) 
-0.2(6) 


-3(1) 
-4(2) 


0.11(2) 

0.12(2) 


P' 


P' 


P' 


P' 


P' 




P' 

as 


P' 

r ah 


0.9(1) 
0.9(1) 


0.4(1) 
0.4(1) 


0.9(5) 
0.9(6) 


-2(1) 
-3(2) 


0.4(4) 
0.3(5) 


-2(1) 
-1(2) 


-1(2) 
1(3) 


0.08(5) 
0.06(7) 



Table 5: Values of the SU(2)-HMChPT Gt parameters from Fit I of eqs. (El and tfH) 
(upper table) and from Fit II of eq. If21\) (lower table), as obtained by including (first row) 
or excluding (second row) the data at [3 = 3.8. Quoted errors are statistical plus fitting 
errors. 



where the coefficients D[ and D' 2 are expanded as linear functions of the strange quark 
mass as in eq. ffT9]) . 

We find that the HMChPT parameter g c cannot be determined from the fit, which is 
almost insensitive to it. It is thus constrained to the experimental value g c = 0.61(7) [2[ 
4"6"] which is in good agreement with a recent unquenched lattice determination, g c = 
0.71(7) [37]. 

The values of the coefficients Di and D\ as resulting from the fits of eqs. ({TBI) and (|20!) are 
collected in Table [5j Using these results and inserting in eqs. ([TBI) and (1201) the experimental 
values of the relevant meson masses we obtain for fo, fo s and /d s //d the results given 
in Table [6] labelled as SU(2)-HMChPT, Fit I. As in the light meson case, we show in the 
Tables the results obtained by including or excluding the data at (3 = 3.8. The values of 
the fit parameters are consistent in the two cases and the results for the decay constants 
are essentially equal. 

The alternative approach we considered to determine the D and D s meson decay con- 
stants is based on the study of the ratios Ri and R 2 defined in eq. (fTBj) . The SU(2)- 
HMChPT predictions for these ratios are easily obtained by dividing the expressions (|18p 
and (JUD for fay/raj^ and R by the SU(2)-ChPT predictions © and © for f w and f K , 
corrected for FSE as in eq. (JBJ). The resulting expressions read 



Ri 



Pi_ 

1/2 

r 



1 + -Xu In 6, + P 2 £« + (P a + P as £ ss ) -j + P ah a 2 m 2 Ds 



+ 



3/2 ' 

Tq m Ds 



(21) 



R 2 = P[ 



1 + ( 2(1 + 3g 2 c ) - ~) & In& + P&i + (P a + PLU ^ + PL a 2 m 2 Ds 



+ 



P' 



r m Ds 



where the coefficients Pi, Pi, P[ and P' 2 are then expressed as linear functions of the strange 
quark mass as in eq. (fT9"|) . 
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} D [MeVJ 


} Ds [MeV 


£ IP 

3D J 3D 


X /do J 


SU(2)-HMChPT 
Fit I 


incl. (3 = 3.8 
excl. = 3.8 


195(6) 
194(8) 


242(3) 
239(5) 


1.24(3) 
1.23(4) 


{93/136; 61/136} 
{25/96; 49/96} 


SU(2)-HMChPT 
Fit II 


incl. P = 3.8 
excl. = 3.8 


199(6) 
195(8) 


246(3) 
243(5) 


1.24(3) 
1.24(5) 


{146/136; 69/136} 
{24/96; 39/96} 


SU(3)-HMChPT 


incl. P = 3.8 


197(6) 


239(3) 


1.22(2) 


371/179 



Table 6: Values of f D , f Ds and f D Jf D as obtained from the SU(2)-HMChPT Fits I and II 
by including or excluding data at P — 3.8. We also show in the last row the results obtained 
by fitting both fu s Jmo, and fDy/^D with their common SU(3 )-HMChPT functional form. 
For each fit, the chi-squared per degree of freedom is given in the last column. For fits 
based on SU(2)-HMChPT, two values of x 2 /dof are displayed, corresponding to fjj a and 
R for Fit I or Ri and R 2 for Fit II. Quoted errors are statistical plus fitting errors. 

The values of the coefficients Pi and P[ are collected in Tableland the results for fo, 
fr) 3 and JdJId are compared to those obtained from Fit I in Table [61 The two fits yield 
results that are in good agreement and with very similar uncertainties. 

In fig. H]we show the dependence on the pion mass square of the four quantities studied 
in Fits I and II. We observe that fo s -s/ttid s (top-left) has a very mild dependence on 
m^, in agreement with the SU(2)-HMChPT prediction of eq. ([TBI) according to which 
chiral logarithms are absent for this quantity at NLO. The logarithmic dependence in 
fu a yjrri£)j fx (top-right), thus comes only from the chiral logarithms predicted by SU(2)- 
ChPT for the kaon decay constant, see eq. ([3]). The lattice results for the double ratio 
Wd s \J m Dsl Ik) I (/dv^d/Ztt) (bottom-right) are almost independent of the light quark 
mass. This is not unexpected, since the chiral logarithms largely cancel in the ratio. We 
also note that in the ratios R and R 2 (bottom plots) , where fjj a ^Jrn Da is divided by fo y/f^D, 
discretization effects turn out to be negligible, smaller than the statistical uncertainties. 

As done for the light mesons decay constants, as an alternative to the SU(2)-HMChPT 
approach, we have also tried to fit both the D and the D s decay constants using SU(3)- 
HMChPT. The corresponding fitting formula is given by eq. f[T7j) with the addition of 
discretization terms. As for the SU(2) case, we have tried two different fits, in which 
foyJ^D and /d s a/"757 are either divided or not by the light decay constants /„. and fx- 
At variance with our results for the light meson sector, we find that the quality of the 
SU(3) fits, with x 2 /dof ^ 2, is worse than in the SU(2) case. For illustration, we show the 
results obtained from the SU(3)-HMChPT analysis in the last line of Table [6] (for the case 
in which the light decay constants are not introduced in the ratios) and in fig. [H Even 
though these results are consistent with those obtained from SU(2)-HMChPT, given the 
poor quality of the SU(3) fits, they are not considered in deriving the final results. 

4.2 Results for f D , f Ds and f Ds /fD 

The results presented in Table [6] show that the SU(2)-HMChPT analyses based on Fits 
I and II lead to determinations of fc, fo a and JbJ fo that are in very good agreement, 
with very similar statistical uncertainties. We choose to average these results and to quote 
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Figure 4: From top-left to bottom-right: lattice results for fo a ■JmrT s , R\ = fD s ^/^D a /fK, 
R = fD s ^mr s /{fD^/m~D) and R 2 = {f Ds ^/rn^/ ' f K ) /{fDy/rnB/ ' U) as a function of the 
pion mass square = m,ps(fJ>i, fa, Hi) 2 , in units of r . We display data with /i s and \i c 
fixed to the simulated masses that correspond to reference strange and charmed meson 
masses romps{fii, }J> a , fJ> s ) = 1-63 and romps{Hh A 4 *) A*c) — 4.41. The SU(2)- (SU(3)-) ChPT 
extrapolation to the physical pion mass is represented at fixed lattice spacing by the dashed 
(dotted) curves, and in the continuum limit by the solid (dashed- dotted) curve. The 
physical results, illustrated by diamonds in the plots, are obtained from SU(2)-ChPT after 
interpolating to the physical strange and charm quark masses. 



their deviation from the average as an additional systematic uncertainty due to the chiral 
extrapolation. 

As in the light meson case, we estimate the uncertainty due to residual discretization 
effects by assigning an error of 50% to the extrapolation from our finest lattice at f3 — 4.05 
to the continuum limit. In the former case we obtain f D ~ Am = 208 MeV and f D ~ Am = 
257 MeV, that are ~ 5% above the continuum limit estimates. Note that this effect is 
larger than the naive estimate of leading discretization effects as being of 0(a s a 2 fi 2 ), which 
follows from the observation that 0((afi c ) n ) effects have been corrected at tree level in the 
definition of the decay constants. This finding clearly illustrates the importance, for lattice 
studies of heavy quarks, of evaluating discretization effects with simulations performed at 
several values of the lattice spacing, rather than on the basis of simple order of magnitude 
estimates. We also find that discretization effects largely cancel in the ratio of the decay 
constants, and we obtain (/^//d)^" 4 ' 05 = 1-23 from Fit II, that is only 0.8% below its 
continuum limit estimate. The same difference is even smaller in the case of Fit I. 
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Figure 5: Lattice QCD determinations of the D-mesons decay constants fo s (top) and fo 
(bottom) obtained from simulations with Nf = 2 [49} I50f and Nf = 2 + 1 \4l\ I51f dynamical 
fermions. A star in the legend denotes preliminary results. The lattice results for fu s are also 
compared with the PDG 2008 experimental average t '2| and with the recent improved measurement 
by CLEO J2f. For fo we compare with the CLEO determination [48]. 



Our final results for the D and D s decay constants and the ratio Jds/Jd are then 

f D = 197(6) (2) (6) MeV , fjj a = 244(3) (2) (7) MeV , f D J f D = 1.24(3) (0)(1) , (22) 

where the errors come from statistics plus fitting, chiral extrapolation and discretization 
effects, respectively. By combining all these uncertainties in quadrature we finally obtain 

/d = 197(9) MeV , fo s = 244(8) MeV , f D Jf D = 1.24(3) . (23) 

The result obtained for fjj is in very good agreement with the CLEO measurement, f^ p ' = 
205.8(8.5) (2.5) MeV [48], and with other N f = 2 03 ED] and JV> = 2 + 1 gD EI] lattice 
calculations, as shown in fig. El Even more interesting is the comparison shown in the same 
figure between our result for fu 3 , other lattice results and the experimental measurements. 
The PDG 2008 average was f^' = 273(10) MeV 0, higher than the values indicated 
by lattice calculations, for which a possible explanation as an effect of new physics was 
given in refs. [211 122]. Recently, however, CLEO has performed with higher statistics 
an improved measurement of the branching ratio Br(D+ — > t + v — > e^vvv) [52] which, 
combined with their measurements of Br(Df — > and Br(D+ — > t + v — > 7i + uu), gives 
fo s = 259.5(6.6) (3.1) MeV [3]. This latter determination, being in better agreement with 
our and other lattice results, weakens the possibility of a new physics effect in leptonic 
D s decays. In ref. [3], also an improved determination of the ratio of D s and D decay 
constants is provided by CLEO, /d s //d = l-26(6)(2), which is in good agreement with our 
result in eq. (T2"3"j) . 
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